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' We focus on the problem of an impurity-free billiard with a random position-dependent boundary 

, coupling to the environment. The response functions of such an open system can be obtained non- 

■ perturbatively from a supersymmetric generating functional. The derivation of this functional is 

based on averaging over the escape rates and results in a non-linear ballistic cr-model, characterized 
by system-specific parameters. Particular emphasis is placed on the "whispering gallery modes" as 
the origin of surface diffusion modes in the limit of large dimensionless conductance. 
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Wave billiards with smooth boundary walls find a wide variety of applications in condensed matter physics P|, 
^2 , microwave and acoustic chaology 0. For chaotic transport in nano-size two-dimensional electronic systems such as 
quantum dots, wires, junctions, and corrals, billiards serve as convenient theoretical models for a confined electronic 
do ' gas. Chaos in a cavity might be attributed to the presence of impurities or to the properties of the boundary itself. 
^ . A review of closed billiard dynamics can be found in Jj3j. Convex hard- wall cavities open to the environment are also 
well studied within the semiclassical approximation 
. An intermediate situation arises when one allows for a small but finite coupling to the outside world all around 
^ ' the boundary. Furthermore, if this coupHng is described by a random set of_coupling coefficients to a large number 
^ [ of ideal leads, the resulting system acquires a natural statistical description Q. Indeed, in studies of nanostructures 
' . > one typically aims at the statistical properties of various response functions, computed through the cavity Green's 
^ ' function. In experiments the statistics are obtained via spectral averaging or from an ensemble of different boundary 
Q . configurations, corresponding to the same spectral density. In both regular and chaotic cavities, the most well studied 
Q ' statistical ensemble is the one defined over the different configurations of the internal impurities 0, Q, Bufl] • 

' , The notion of a statistical ensemble does not come naturally for the analytical study of a biUiard with no internal 
' scattering potential, no magnetic impurities or random magnetic field. The scattering, which in this case takes place 
^ , exclusively at the boundary, is the source of stochasticity. In the simplest case the scattering is specular, but it can 
\^ ■ also be diffusive, as in Furthermore, the scatte ring may be accompanied by the escape of scattered waves into 
OO ! the exterior fF|, as, for example, in quantum dots or quantum corrals In this paper we present a more 

. detailed description of the method proposed in 0. The model system we study is a biUiard weakly attached to the 
large number of ideal infinite leads For simplicity and clarity of presentation we consider one open channel in each 
waveguide, with a coupling coefficient which is Gaussian distributed around a small, but non-zero mean. The mean 
and the width of the coupHng coefficient distribution are assumed to vary smoothly with the position of the lead on 
the perimeter. 

' For this dot, we construct the supersymmetric generating functional 0, 01 and performing the average over the 
I ensemble of realizations of coupling coefficients, we obtain a "surface" non- linear cr-model (NLcrM). Its "diffusion" 
I ' modes are confined to the boundary of the dot. In convex nearly closed biUiards these correspond to the so called 
whispering gallery modes (WGM), which are high angular momentum modes corresponding to the classical trajectories 
C . running alongside the billiard walls. The WGM are exponentially less likely to escape compared to the modes with 
incidence directions close to the lead normals, as can be inferred, for example, from Thus, the response functions 
at very large times are expected to be dominated by these modes, and hence can be non-perturbatively calculated 
from our supersymmetric functional. 

Some of the technical tools we employ, such as the non-Hermitian effective Hamiltonian have been frequently used 
along with random matrix theory (RMT) 0,0,0]. The use of RMT impHes a restriction to the universal regime. 
] Here we pursue a description of the nonuniversal regime of the billiard dynamics, thus covering a much broader 
energy range than is possible with RMT. We are also interested in spatially dependent characteristics of the system, 
inaccessible via RMT. 

The assumption of only one open channel is not essential. The generalization of our approach to an arbitrary 
number of open channels is straightforward. The choice of random Gaussian distribution of the channel strengths is 
in fact quite realistic; the actual values of couplings (and the number of open channels) are not known in general [T^ . 
It is usually argued, at least in situations where the RMT is applicable, that quantitative results do not depend on 
the coupling being random or constant 0|. For our development, however, the randomness is essential as we briefiy 
explained above. 
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We note that our model applies, mutatis mutandis, to a system with local or extended sources of damping |l7lll8l|. 
Besides quantum dots there are other mesoscopic billiards which fall into this category, such as quantum and 
optical corrals optical resonant cavities and the artificial atoms proposed in Here also we expect the 

WGM to play an important role in the long-time Green's functions. The dependence of escape rate on the angular 
momentum, important in both optical and electronic systems, is also incorporated in our model. 

Finally, weak interactions in nanostructures with large dimensionless conductance creates an additional motivation 
for our study. For a ballistic cavity, in which electronic interactions are of interest, while the shape is not, the convex 
smooth- wall billiard represents a sufficient starting point. It might be possible to use the results of our present 
analysis in the large- approach of to the interacting ballistic case. We also note in passing that the fundamental 
problem of constructing a a model for a closed, impurity free ballistic billiard has proven to be technically challenging 
[lO, 23, 24, 25]. Our work shows that a random coupling to the environment acts as a natural regularizer and allows 
us to circumvent many of the technical difficulties. 

The plan of this paper is as follows. Section II reviews the procedure for "integrating out" the leads to arrive at a 
non-Hermitian Hamiltonian which describes electronic scattering in a quantum dot. In Section III, we construct the 
supersymmetric generating functional for the correlation functions and carry out the ensemble average, introducing 
a supermatrix field for each lead to decouple supervector fields. The continuous model resulting in the limit of large 
leads number is analyzed in Section IV. We demonstrate how a Born-Oppenheimer-like approximation j2^ can be 
employed in the solution of the effective problem, which results from the saddle point evaluation of superintegral. 
This is done in Section V. We complete the derivation of the surface diffusion NLctM in Section VI and present our 
summary in Section VII. 



II. EFFECTIVE NON-HERMITIAN HAMILTONIAN 



It is common in mesoscopic physics to study open systems using their closed counterparts and vice versa. The 
eigenfunction structure of closed dots affects the transport properties of the open dots, while the leads attached to the 
dot can change the nature of the dynamics from regular to chaotic [l^- In the RMT context, the scattering approach 
[T^ to mesoscopic billiards relies on the decomposition of the system into internal and external subparts. The internal 
is described by N bound states, while the external, by M channel states propagating along the hard wall leads to 
infinity. To set up the analysis one introduces (i) the N x N internal Hamiltonian Hin, and (ii) M x M ^'-matrix, 
which relates incoming and outgoing wave amplitudes in the asymptotic region: 



J,. _ X.., /.incoming , r<(ji) /Outgoing 



where z, j specify leads, and m, n specify channels, (iii) express ^-matrix in terms of Hin and N xM matrix W , which 
couples the subparts S ^ I — 2TTiW'' {E — Heff )~ W, with Heff — Hin — ittWW^ . This relation of S'-matrix to the 
effective Green's function {E — Heff) ^ is the building block of the Hamiltonian approach to the system's statistics 
in the universal (RMT) regime, i.e., the regime independent of the details of the underlying classical dynamics [T^ . 
The non-Hermitian random matrices have been the subject of numerous works on quantum chaotic scattering (see, 
for example, and references therein). 

Similar steps lead to the generalization of the effective non-Hermitian Hamiltonian in the nonuniversal regime, as 
was done in a thick wire i n 12^ and for a disordered quantum dot, in which electrons, confined by a hard wall potential 
can escape into leads, in (2g |. 

In full analogy with the RMT problem one can take the sam e ap proach to our problem involving open ballistic dot. 
For completeness, we reproduce the main steps carried out in j22|, noting that we deal with a two-dimensional clean 
electronic system. 

In order to study the Green's functions G 

{E±ie- Ho) G^'^ (r, r') = 6 {r - r') , (1) 
corresponding to the dot-plus-leads Hamiltonian 



e 



c 
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Ho= P--A . (2) 



we introduce the crossectional surfaces C„ perpendicular to the walls of the leads close to the place of their attachment 
to the dot. Then we reformulate the problem by specifying the boundary conditions on these surfaces, thus eliminating 
the leads at the expense of modifying the Hamiltonian of the system. First we introduce auxiliary retarded and 
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advanced Green's functions G ' (r, r') for each lead, satisfying Eq. inside the region of the lead bounded by the 
section C„. The boundary condition at C„ is the only difference between G ' (r, r') and G^''*(r, r'). The former 
vanishes at this boundary. Considering G^'"^ (r,r')j such that r and r' belong to the regions on the different sides 
from the section C„ we make use of the following identity (current conservation) 

G^ (r, r")) * v,G^ (r, r') + G" (r, r') (v,G^ (r, r")) *] = 2iS (r - r") (r, r') , (3) 

where = —iiVr + eA/c) /to. To arrive at the current conservation relation (Eq. ijSl) we combine definitions of 

both G^ and G (Eq. (QJ), by multiplying them with G^ (r,r') and subtracting from each other. Then, integration 
with respect to r-coordinates taken over the volume of the lead results in 

G« (r", r')^-U (x„ . V.) (g^ (r, r")) * G« (r, r') dr, (4) 

^ 

Upon application of boundary condition on G (r, r"). Here is the unit normal to Cn. Next, we apply the velocity 
operator Vr" to both sides of Eq. Q and pick r" = y on the crossection Cn to get 

(x„.v,)G«(y,r')= / i3„ (y, r) G^ (r, r') dr, (5) 



where 



2 

(y, r) = - (x„ • Vj^) (x„ • Vr) G (y,r). 



and we used (G (r,r")j =G (r",r). Note, that the functions -B„ (y, r) are completely determined by the properties 

of the leads, which we assumed to be ideal. Thus, the "dot-plus-leads" system governed by Eq. (QJ is now reduced to 
the dot only, governed by the same Eq. (QJ and boundary conditions specified by Eq. ij^l. We also note, that it is 
real part of a function Bn (y, r) that is related to the fiux carried into the leads by the electronic states. To simplify 
the boundary conditions we modify the Hamiltonian of the system even further according to 

. N 

H ^HoT^Y.^^^^^' (6) 

ji=i 

where SCn is a surface (5-function, defined via J (5G„^' (r) dr — J^, (r) dr, and -Bn^* (r) — 5R_B„ (r, r') (r') dr' , 
with (r,r') = jnVn'fin (r) (fin (r')- To compacticify the formulas of the next sections we also assume that each of 
the N leads contains at most one open channel. Thus, we end up looking for the Green's functions of the effective 
problem 

{E-H)G''-^{r,r')^6{r~r'), (7) 
associated with the Neumann boundary conditions 

(x„.v,)G«'^(r,r')|c„+o = 0' (8) 

where the derivatives are taken from the side of the lead (which is indicated by Gn + 0). The second term in 
Hamiltonian (Eq. describes the finite escape probability for the electrons coUiding with the boundary. Now we 
are ready to proceed with the construction of the cr-model along the conventional lines developed for closed systems 
0. Note, that boundary condition given by Eq. ^ ensures that system is effectively closed, i.e. component of the 
current, which is normal to boundary, vanishes. 



III. SUPERSYMMETRIC GENERATING FUNCTIONAL 



Next we construct the generating functional for both retarded and advanced Green's functions for Gaussian dis- 
tributed dimensionless coupHng coefficients 7„. These coefficients are related to sticking probabilities and transmission 
coefficients, frequently used within the Hamiltonian approach to chaotic scattering to compute statistical distri- 
butions of the resonance widths, delay times and related characteristics. A detailed discussion of the physical meaning 
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of these coefficients can be found in 0, ItR and in references therein. In order to perform non-perturbative 

calculations of these averages in the non-universal regime we first construct the supersymmetric functional Z [flLl3ll|. 
Any correlator of the dot Green's functions can be later obtained from it by differentiation with respect to sources J 
as long as we know the average Z [J] over the 7-ensemble. 

We can assume the magnetic field (Eq. to be vanishingly small and remove it from consideration. Its 

role is reduced to breaking the time reversal symmetry, and justifying the use 4-component supervectors 5* (r) = 
{S'l (r) , xi (r) , S2 (r) , X2 (r)} in the supersymmetric functional. The generalization of our formalism to other sym- 
metry classes is quite straightforward and will require doubhng the space It is convenient to express the coupHng 
coefficients as a sum of constant and stochastic parts: 7n = 7„ + Jn- For the statistics of 7„ we assume that (7„) = 0, 
i^nlm) = x'^iSnm', all higher moments factorize into second moments. We indicate averaging over random couplings 
to the leads by the shorthand notation {...)-. Then we "eliminate the leads" passing to the Hamiltonian given 
by Eq. ®. 

In terms of supervectors ^' (r) , and supermatrices L — diag {1,1 — 1,1}, A = diag {1,1 — 1,-1} the generating 
functional Z [J] is written down as follows 

{Z [J])^ = J d^*d^e~^^^^ {'-'^'^"X ' (9) 

r ^ 1 ^ r 

£ [VI/] = z / vl/t (r) HjL ^ir)dr+-J2 Vn% / (Vn) {Vn) iv'n) AL * (y;) , 



with 



Hj = ( -^-E] h + ieA + J, 
2m 



where e is infinitesimally small, Lpn (y) = y/2/dn sin {ny/dn) (for hard- wall lead of width d„). Here stands for a 
double integration over ?/„ and y'„ , the transverse coordinates along the crossection C„ . The exact form of the source 
supermatrix 

J ^ diag { Ji (r) , J2 (r) , Ji (r) , J2 (r)} , 

is dictated by the choice of the physical quantity we eventually wish to calculate. 
Averaging over 7„ produces: 
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= y + lJ2 { £ (y") (y") * (2^")} +■■■/ ■ (10) 

This step towards constructing supersymmetric NLcrM for our system is no more difficult than in the problems 
requiring averaging over impurities or random magnetic field. As we shall see below it results in the coupling of the 
effective field with the source of boundary transmission. 



IV. EFFECTIVE PROBLEM FOR A SMOOTH CONVEX BOUNDARY 

Next we carry out the Hubbard-Stratonovich transformation to decouple the "interaction" terms (Eq. ifTIl l which 
entered the action of Eq. (pj after averaging. The procedure, involving supermatrix fields Qn, is explained in Appendix 
A and leads to 

{Z[J]),^ f Y[DQ^ I dvI/*d*e-^[*lexpE(^^ / ^Hy^) L 

n n=l ^ ^ "'C„ 

X ^„ (y„) (y;) * (y;j / g„ {yl, y'^) ^„ (y^;) {y'^) 



2 



2 
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Next, we make further simplifications by setting dn = d and distributing the large number of narrow leads dn = d P 
(where P is perimeter of the dot), close to each other all around the boundary. In this limit we can approximately 
write (fin iVn) iVn) — 2/d- These two model assumptions are not essential for subsequent analysis, although they 
do make it more transparent. The distribution of the leads coupling coefficients can be chosen piece- wise continuous, 
while widths can be given as a set of input parameters to the problem at hand together with 7„ , Xn and w„ . 

Then we make use of the mean value theorem for each of the integrals, provided all the integrands are smooth 
functions of y„. The resulting sum in the exponent of Eq. Ijllll can be combined into a single integral over the 
arclength s, ranging from to P, which runs along the boundary, defined in polar coordinates as r = i? (s), producing 

{Z [J])^= j DQ{s) j d^* {r,0) {r,9) e"^!*! 

X exp 1 2md [ x (s) {R (s) , s) LQ (s) {R (s) ,s)ds- 2m^d [ (s) ds \ , 



where x (s) = x (s) v (s) and both band velocity v (s), mean 7 (s) and r.m.s. x (s) of the coupling strength are now 
smooth functions of arclength. The corresponding changes in the term C [^] are made accordingly. Upon carrying out 
the ^'-integration, we arrive at the following representation of the generating functional in polar coordinates r = (r, 9) 

{Z[J])^ = j L>QexpStr ^2m^d j dvdv' {s) 5p5 [v - v') 

- In { (-.Ho + ^^^A^p) 5 (r - r') - (s) S,S (r - r')} 

- In ih + eAG (r, r') + G (r, r') J (r))] , (12) 
where (5p is a perimeter delta-function and the supermatrix "Green's function" G is determined from 

|ho - *2mi is) d (q (s) - ^-2ifI_A^ Sp^G (r, r') = iS {r - r') , (13) 

and Ho ~ TCj {J — 0) . Thus we reduced the generating functional for our system to the integral over the supermatrix 
superfield with boundary support. We have 

{Z[J])^ = jDQe^i'^^+^-'i^K (14) 

with the free energy 

F [Q] = Str j drdv' ^-2m'^dQ [sf Sp6 (r - r') + In -iG^^ (r, r')} , 

and symmetry breaking terms 

Fj[G] = -StrJ \n{h + tKG (r, r') + 5 (r, r') J (r)) drdr'. 

In order to reduce our generating functional (Eq. Ijl4|l l to a NLtrM we employ the saddle point condition, which in 
our case reads 

Qsp (s) = (i? is) ,R{s), s, s, Qsp (s)) . (15) 

We assume the saddle point solution to be diagonal: Qsp (s) = Qo (s) A. To proceed with the analysis of fiuctuations, 
one needs to determine both the Qq (s) and the diagonal Green's function supermatrix Gsp (r, r')- Thus, by combining 
Eqs. lfT3jl and itTfill with the assumption about saddle point structure we mapped the original problem with random 
boundary condition onto an effective problem specified by the differential equation: 

1/9 9N^/ 9N iS(r-r')5(e~e') 

2m r 
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where = —2mE, with associated boundary conditions 

|-g (r, r\9, 0', |,„ = (r, /, 0, 0', |,_ , (17) 

^g{r,r',e,9',n')\^.,=0, (18) 

where / (s) = m^a; (s) Qo (s) (s) d, , 5+ are the inner and the outer surfaces of the dot, and Qo (s) — Qq (s) — 
7(s)/(2x(s) md). 

V. BORN-OPPENHEIMER-LIKE APPROXIMATION 

To construct the Greens function of Eq. Ijl6|l we employ the technique of j^, which was also used in Q for a 
circular billiard. In the latter case the corresponding Green's function reads 

oo 

gcircieir,r\e,e\^) = — V /„(«;r<){a„/„(^r>)+if„(«;r>)}e™(^-^'), (19) 

n— — cxj 

with r> (r<) is a maximum (minimum) of r and r', /„ and Kn are modified Bessel functions respectively, and the 
coefficients a„ are chosen to ensure the boundary condition. The summation in Eq. Ijl9|l is replaced with integration, 
while Bessel functions are replaced with their uniform approximations 0|. 

For a generic convex smooth-wall billiard a similar expansion is possible with the help of the recently proposed 
Born-Oppenheimer-like approximation j2^. The method used in assumes that the angular variation of x{s), 
f (s), etc. is slow and enables one to determine certain classes of eigenstates quite accurately. Guided by the example 
of circular billiard, and the fact that WGM have the longest life-time, we focus our attention on the limit of large 
angular momentum. 

The Greens function (Eq. (QHl) was build from two radial solutions of Eq. Ijl6|l and angular harmonics. In case of 
a generic billiard, the angular "part" of the eigenstate is obtained from WKB solution of slow- variable equation (2^ . 



r f F (s) 
V'(s) = exp|jAy 'r(I)^^ 

while the radial "part" is given by the Hnear combination of /((s) (rer^) and if;(s) (urs)- Here angular momentum 
"index" is parametrized as I (s) — XF (s), is a local radial coordinate, computed from the local center of curvature 

[2^ 1. Eigenvalues A are determined from arclength quantization condition (-0 (0) = tJj {P)): A„ = 27m/ (^^^^P- Here 

and below we use "overbar" to denote the perimeter average : R (s) ds/P, f ^ f (s) ds/P, etc. Note that 

for dimensionless quantities we will not use overall factor l/P. The explicit form of F {s) can be found from the 
radial quantization condition Below we consider F (s) to be known. Thus, the Green's function for the effective 
problem given by Eqs. Hl(ill7|l can be approximated as 

n— — oo 

X exp \iXn f ^^dA , (20) 



with 



a„ (s) 



Ris) 

-if (s) Kx„F(s) (Ki? (s)) + kR (s) K'^^p(,) {nR (s)) 

(s) h„F[s) i^^R (s)) - (s) /;„j.(,) («-R (s)) 

the boundary, but it suits our purposes. Fi 
We have, at the billiard boundary 

imR{s) I\„F{s) {i^Pjs)) 



This expression does not hold far from the boundary, but it suits our purposes. For example, it can be used to 
determine saddle point solution Qo (s)- We have, at the billiard boundary 

g{R {s),R{s),s,s',K^) 



X exp i iXn I ds 



~^oo ^A„F(s) {kR (s)) - kR (s) /^^^(^) {kR (s)) 

RisY 
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Then, Qo (s) is determined by the stationary point condition (Eq. dJ) 

.x{s)Ris) ^ /A„f(,)M(s)) 



z- 



E 



Dropping the imaginary part of Qo (s), we set g (s) = kR{s), and proceed as follows. We evaluate the sum over n 
in Eq. lf2T]l asymptotically in the limit: g (s) ^ 1, / (s) /^(s) ~ 1 by replacing it with integral over corresponding 
continuous variable. Following the technique, described in Ref. j^^l) and illustrated in Appendix B for a more involved 
calculation, we use uniform approximation for the Bessel function and its derivative [231) carry out the integration 
and rewrite Eq. lf2T]l as 



x{s)R{s){F{s)/R{s)) ~ 7(s) 
- Qo (s) 



oj;,/ N / f / \2 / n2 2x(s)md 
(s) Y / (s) - g (s) ^ ^ 



after the substitution g (s) —ig{s) {k —ik), to the leading order in l/g{s). This equation can be now solved 
numerically yielding / (s) (or Qo (*)) for any prescribed set of parameters of the problem. 



VI. BALLISTIC NONLINEAR cr-MODEL 



Having determined the relation between the Green's function G and saddle-point value of superfield Q, we come 
back to further analyze the generating functional specified by Eq. II14II . As a final step in derivation of NLcrM we 
have to expand the action F [Q] up to quadratic order around the extremum and include first variation of symmetry 
breaking terms in Fj [Q]. Observing that F[Qsp] = we turn to the fiuctuations of Q, which decompose into a 
transverse piece (SQ^*^ (along the saddle-point manifold 0) and a longitudinal piece JQ^'^ (orthogonal to the saddle 
point manifold). We focus on the transverse part of the action, the part anticommuting with the A- like saddle-point 
solution. Our goal is to demonstrate that in absence of symmetry breaking terms transverse fiuctuations are massless 
(Goldstone) modes of the theory. The purely transverse terms are given by P*] 

Ft[6Q]^ -Am^dStr (^S Q^''^ (s) j ds + {mdf Str / G (Qo (s)) G {-Qo (s)) 

xx{s)x is') 5Q(*) (s) ^Q(*) (s') dsds'. 

In view of the developments of the previous Section, we can expand deviations of Q in approximate angular eigenstates: 
SQ^*'^ (s) = Z^i^-oo Qf^ 6xp {iXi J'' dsF (s/i? (s))}, and, setting G (Qo) = G, and G {-Qo) = ^ we obtain 

Ft [6Q] = -4m^d Str ^ <50[')(5gWe^(^'+^") ^"^ds + {md)^ Str ^ SQ^^SQ^^^ 

l,n •'^ l,n 
nP nP _ , 

X ds ds'GGx{s)x{s')e'^'-f^ (22) 
Jo Jo 

Our intermediate goal now is to get a Ward-like identity, which would allow us to evaluate massive part of the 
transverse action. At this point we consider Eq. lfT3jl for both Green's functions. 



and multiply Eq. ll23ll by Gh{s,s') and Eq. 1I2H1 with Gh{s,s'). We are free to chose function h{s,s') to behave 
as we like inside the billiard, because only the value of this function at the boundary matters for manipulations in 
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Eq. H22|l . In particular we may require that h slowly decay to zero in the radial variable as we move away from the 
boundary. Since the radial rate of variation is approximately l/R, while the angular rate is 1/Af (where A_f is a 
Fermi wavelength), we neglect all radial derivatives in the subsequent analysis, which is correct to leading order in 
R/Xf- Next, we subtract one equation from the other, integrate the resulting expression over the area of the dot D 
(including the boundary: < r < i?+ (s)), and set the remaining free radial coordinate to R{s) to get 

J {hgv -vg -hgv -vg^dT dsf {s)h{s,s') 

xg (i? (s) , R (s') , s, s')g (i? (s) , R is') , s, s') + 2imh (s', s') 
X (g (i? is') , R is') , s', s') ~g (i? is') , R is') , s', s')) . (25) 

Note, that according to the saddle-point condition (Eq. Ijl5ll ') 

g iR is) ,Ris),s,s)-g iR is) , i? (s) , s, s) = -2m£^, (26) 

X [s) 

which might be used in the last Hne in Eq. ll25|l . Furthermore, one can transform the left-hand side of the Eq. lf25|l 
by making use of Gauss's theorem for the domain D (with boundary dD) 

J (hg\i ■ vg - hg\/ -vg^dr^ J h (gvg - g\/g^ ■ da 

- y" (V (/le) ■ vg - V ihg) ■ vg) dr. (27) 

Application of the boundary conditions (Eq. Ijl8|l l to Eq. Il27|l makes the surface integral vanish. It remains to 
evaluate the area integral on the right-hand side of Eq. ll27ll . which we now substitute into Eq. (l25)l and obtain, 
neglecting the terms with radial derivatives of h 

= 2i [ dshis,s')gg, (28) 
Jo 

where we used Eq. ll26)l . The boundary value of function h is completely at our discretion. We select 
h is, s') = -^ imd)' Y: ^MlMe^^' ^ 0}'^+^^" /=' ^'^^qW^qW. 

Then we integrate both sides of Eq. Ij28|l with respect to s' along the boundary of the billiard and substitute the 
result together with h in each term of the sum in the second term of Ft [SQ] (Eq. 12211 1. 

We observe that the massive contribution coming from variation of Q (s)^ (the first term on the right-hand side 
of Eq. 1I22I1 I will be canceled by the second term on the left hand side of Eq. lf28jl . The surviving linear in JQ^*' 
term proportional to 7 (s) (cf. 0) will be accounted later, together with the symmetry breaking terms. In addition 
to that, we need to demonstrate that massive longitudinal modes (SQ*^'-* are decoupled from the essentially massless 
transverse modes. Then, the longitudinal fiuctuations around saddle-point can be all integrated out, producing unity 
due to supersymmetry. 

To complete the derivation of NLctM we still have to analyze the remaining part of the transverse action 



Ft [6Q] - Str V [ ds' [ ds [ 



^^'^ dh is,s') ^ , , /- d ^ ^ d ~\ dr 

^ ' ^R^is){g—g-g—g]— (29) 



ds \ ds ds' 



We perform this step using the asymptotic technique, we employed in previous Section. After a straightforward, but 
lengthy procedure, which we allocate to the Appendix B we obtain 
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"4^1 f'^P~9^ 



One can use similar set of manipulations for the purely longitudinal and proportional to SQ^'''' SQ'^^^ + ^Q'^'-'cJQ'^'^ parts 
of the action, which consist of the fluctuations commuting with the saddle-point solution, and therefore contain Q'^ in 
the contribution to second variation coming from the logarithmic term in F[Q] (cf. Eq. lf22|l ). These manipulations 
enable us to show that to the considered order in I/5 the longitudinal and transverse modes are decoupled (the 
coefficient in front of the JQ^'-'i^Q^*) + SQ'-^^SQ'-^^ is small compared to the mass of the longitudinal modes). 

To finish the construction of the nonhnear cr-model, we integrate out the longitudinal modes, set SQ^^^ — Q, and 
expand the symmetry breaking terms Fj [Q {Q}] to the lowest order in J and e to get 



{Z[J])^^ J DQe 



F[Q] 



with the free energy given by 

F [Q] = Str / drdr' 



where 



Dn 



R dQ (s) 
f ds 



xmd 



'{s)K\5p5{y-v') 



J ds" (eA + J (r)) Q (s") a (r, R (s") s"; r', R (s") . 



(31) 



mdx 



a (r, R is") s";r\ R {s") s") = t^Gsp {R {s") s" , r') (r, R {s") s") 

and we suppressed explicit slow s-dependence everywhere. The free energy given by Eq. l(3T|l is the central result of 
our paper. It displays the surface modes Q (s), which undergo diffusion and drift, and are coupled to the interior of 
the dot by the last term. 

For the calculation of any physical quantity, expressed via Green's function correlators such as, e.g. 
(G^ (r, r') (r, r'))^i it is necessary to have actual parametrization for Q (s). Since basic symmetries of our problem 
are not any different from these of the corresponding Q (r)-field of the diffusive problem 0,0,1^ our supermatrix Q 
can be parametrized as suggested in 



VII. DISCUSSION 



We have constructed a nonperturbative framework to analyze one particular realization of a whole class of nanos- 
tructures: a nearly closed system with ballistic internal dynamics and random losses at the boundary. Our approach 
uses a natural regularizer, which enables us to circumvent the technical difficulties of previous approaches to closed 
baUistic systems [l^, l2Ji U% 12^ ■ We find that the resulting theory can be characterized by diffusive modes confined 
to the boundary and interacting nonlocally with the interior as encapsulated in our main result (Eq. iPIJ)- These 
diffusive modes are identified as WGM, which are exponentially long-lived compared to other trajectories scattering 
off the boundary non-tangentially and, therefore are anticipated to dominate the long-time behavior of response func- 
tions. We expect our formalism to be useful whenever a ballistic nanostructure supports such modes. Note that no 
RMT assumptions have been made, and our approach is applicable in a broad energy range and includes the spatial 
dependence of physical quantities. 

We derive our results by treating a baUistic biUiard as a scattering system. In particular, we use an effective non- 
Hermitian Hamiltonian, which is a common tool in quantum chaotic scattering. However the resulting supersymmetric 
NLcrM is fully intended for the calculations of "internal" characteristics, such as local density of states statistic or 
correlators of eigenfunctions, rather than statistics of refiection or transmission. For that reason, the escape rate is 
chosen to be small and distributed along the boundary. In addition, we chose the coupling part of the non-Hermitian 
Hamiltonian to be fluctuating, thus deflning an ensemble of billiards with different escape rates. These steps allow 
to study ballistic nanostructures without the introduction of any other regularizer 0, 01 • The simplest realization of 
such model is a circular billiard with constant parameters, for which the surface diffusion NLctM was obtained in 
using the saddle-point solution. This biUiard is also the classical example of a system supporting WGM. 

In this paper we have extended the results of Q to show that surface diffusion modes are present in generic, chaotic, 
smooth-wall billiard as well. To create an extension to other chaotic optical and electronic systems, for which WGM 
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are of importance [H El S , one needs to adapt our procedure to the details of the coupling to the environment. 
However, for most mesoscopic structures the wall absorption, which is always present in experiment, can be adequately 
modeled by uniformly distributed leads. Furthermore, our model allows us to add finite number of wide leads if the 
problem under consideration involves transmission properties. Other calculations within a similar framework can be 
performed for systems which have disorder concentrated at the boundary, e.g. quantum corral composed of different 
atoms or artificial atoms proposed in j^J] with the random magnitude of the magnetic flux. 

As for the handling the calculation of correlators of physical quantities one needs to introduce a difference between 
two energies u into the derivation. For the final expression for the free energy (Eq. (EH) this means a replacement 
e ^ e + iuj /2. A good starting point can be the application the perturbative technique described in Ref. [23 to the 
analysis of the 2-point function in a circular cavity. Assuming Q [9) (0 being the polar angle) to fluctuate weakly 
near the origin QqK one can decompose the supermatrix Q {9) as T~^Q {9)T, where T are angle independent 
pseudounitary supermatrices, which are, in turn, parametrized in terms of off-diagonal supermatrices W The 
perturbation scheme is based on the expansion of Q (9) around the saddle-point solution (QqA) in terms of W (see 
the details in jsll). 

It is also natural to inquire about connection of our model to the existing ballistic analogues of the diffusive NLtrMs, 
e.g. the models proposed in 0,E,[23|- Note that the zero-dimensional versions of these theories produce RMT, which 
makes the answer to this ques tion important for studies of truly ergodic systems. The construction involving diffuse 
boundary scattering (Ref. |inl|l is essentially different from ours, because in Ref. the electron loses memory after 
a single boundary coUision, while in our model, the WGM trajectories retain phase coherence until the electron flnally 
leaves the system. Our treatment is complementary (only in conceptual sense) to that of Refs. 0,1131, in which only 
modes of Q inside the ballistic dot appear, and the boundary value of Q is nonfluctuating. The zero-dimensional case 
Q (r) — const inside the dot {Q (9) = const in our case), i.e. the situation in which only lowest mode contributes, 
seems to be the only possible connection between these models and ours. However, another crucial difference between 
the two formulations is worth mentioning. The universal parameter r (mean free time between collisions) in the 
Muzykantskii-Khmelnitskii cr-model (^3|) is angular momentum dependent in our case. In order to address this issue 
of zero-dimensional correspondence between our result and these of Refs. 0,113 one needs to carry out the calculations 
for correlation function of a chaotic system at hand using the NLcrM we derived and compare it to the RMT result 
in the appropriate limit. This is a subject of ongoing work to be presented elsewhere. 

Finally, we hope to consider applications of our NLcrM to the interacting-electron problem in the future. One of 
the possible ways to take the interactions into account in diffusive and baUistic systems with large dimensionless 
conductance, is to use a "Universal Hamiltonian" 01, which was shown to be the renormaHzation group flxed point 
for weak interactions [H ■ We hope to extend our analysis to the interacting ballistic case by using the large- TV 
approach of [2^ . 

We are grateful to the NSF for partial support under DMR-0311761. 



Appendix A: HUBBARD-STRATONOVICH DECOUPLING 

In this Appendix we derive the identity which decouples supervector variables in the action averaged over 7„ . 

For arbitrary function V{u,z)l4, which does not possess any supersymmetric structure and for two coordinate 
dependent supermatrices Q (u, z) and A (w, z) which have the same supersymmetric structure and, therefore commute, 
the following generalization of the decoupling rule can be verified: 

DQe-^^' ^ ^ ""^^^^-'l^^^^^ = Q-Ts:^'^^'' n AVAV ^ ^^^^ 

where we skipped coordinate dependence on y„ and y'^ and corresponding differentials. First of all, we make use 
of the theorem due to Parisi-Sourlas-Efetov-Wegner (see, for example, [2a|), which states, that for any "invariant 
superfunction" the corresponding superintegral is equal to its value at the origin, which in our case is unity. By 
shifting Q ^ Q + i5A, and setting 2a5 — /3 = 0, we get 

-a I QV [ QV + il3 [ QV [ AV >^ - a [ QV [ QV - i {2a6 - 0) 
Jc„ Jc„ Jc„ Jc„ JC„ JC„ 

X J AV J QV^~aJ J QVQV J J AVAV. (A2) 

Then, we carry out the Hubbard-Stratonovich transformation by applying Eq. ijAlll from right to left to each of 
the terms in Eq. ifTTHl (in our case = x'^vfjS and V (yn, y'n) = 'Pn iUn) fn iUn))- We choose a — so that 

if] = mXnVn/2 and get Eq. Illl|l . 
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Appendix B: DIFFUSION TERM 

Approximating the Green functions G [r, R (s') , s, s') and Q (r, R (s') , s, s') according to Eq. Ij20|l and h with its 
boundary value, we carry out r-integration with the help of the identity 



/ 



[z) Ip (z) ^ IL (z) Ip {z) - (z) (z) 
dz — z 



a2 -^2 
and get 

2j\2 rP rP 



d\ ^ f jif , X (s) X {s') g {s) Dp^g {g (s)) Xn 



Ft[SQ]=^(^) E / ds' ds 



P J ,^ Jo Jo /(s)(Ap + A,) 

l,n,p,q \ / \ 1- ^ / 



where Dp^g = 1/ (jp,,/lrp,5 + if {s)g{s)) with 

jp,, = fh^F (g) h,F (5) + q'iLf i~9) iLf (g) , 



tt5p,g = ix^F (g) I'x.F (g) - I'x^F ig) Ikf (g) ■ 

where we omitted s-dependence for brevity in the last two lines; we continue to do that throughout this Appendix. 
Integration over s and s' in the above expression for Ft [5Q] produces infinitesimally small result unless n + p + q = 0, 
I = —n. Anticipating this averages we concentrate on the evaluation of the sum over p and q restricted by these two 
conditions. Later on one can approximate Ft [6Q] by replacing the result of the summation (together with the rest of 
the integrand) with its perimeter averaged value. With this in mind we manipulate J2i npq i^^^o a double sum 

oc 

l,n,p,q l.n— — oo 

|i| \n\ 

The sum over \n\ can be performed asymptotically in the limit: g {s) :s> 1, f (s) /g (s) ~ 1 We convert this sum 
into the integral over the new variable 

271 \n\F{s) 



g{.s)PiF/R) 

use uniform expansion for the Bessel function I\^^^p(^g^ {^X\n\F (s) / fj) f3^ and expand j|n|,|n|+|i|/fD|„|+|;| and 
j|„|j„l_[i|/tti|„|_|;| to the leading order in 1/g . We have 

1 - /i2 - 2fi^ + i^) a: 



2^ (£'t„|,|„|+|;| + L'|„|,|„|_|i|j -— / d^i— , / 77 

|nl Jo g(l + /i2r'(l + M' + |2 



Then, the first integral with respect to fi in Eq. ljBl|l vanishes, while the second one yields 
{jr/'if^g) |(/2 + 252) /■^/2 — g2 — 2,g| after we substitute g(s) —ig{s). To arrive at Eq. ij^fijl we drop some 
of the terms which are negligible in the above mentioned limit. 
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